Shape transformations in rotating ferrofluid drops 
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Floating drops of magnetic fluid can be brought into rotation by applying a rotating magnetic 
field. We report theoretical and experimental results on the transition from a spheroid equilibrium 
shape to non-axissymmetrical three-axes ellipsoids at certain values of the external field strength. 
The transitions are continuous for small values of the magnetic susceptibility and show hysteresis for 
larger ones. In the non-axissymmetric shape the rotational motion of the drop consists of a vortical 
fiow inside the drop combined with a slow rotation of the shape. Nonlinear magnetization laws are 
crucial to obtain quantitative agreement between theory and experiment. 

PACS numbers: 47.20.Hw, 47.55.Dz, 75.50.Mm 



The equilibrium shapes of rotating fluid bodies are of 
importance in various fields of physics as, e.g., astro- 
physics H], nuclear fission [||, plasma and biologi- 
cal physics [Q. The famous controversy between New- 
ton and Cassini on whether the Earth has the form of 
an oblate or prolate rotational ellipsoid started a series 
of ingenious investigations of the equilibrium shapes of 
heavenly bodies including work by Maupertius, MacLau- 
rin, Jacobi, Riemann, Poincare and others. Among the 
suprising results discovered are the possibility of three- 
axes ellipsoids as shown by Jacobi and the emergence of 
pear-shaped configurations found by Poincare. In rotat- 
ing non-neutral plasmas, laser cooled in a Penning trap, 
a novel equilibrium state with non-axissymmetric surface 
has recently been observed Tank-treading elliptical 
membranes in a shear flow have been used as a theoretical 
model for the motion of human red blood cells ||^ . 

Experimental investigations of the stationary shapes 
of rotating bodies in the laboratory usually start with a 
static drop of fluid floating in another, immiscible fluid 
of the same density. The rotational motion of the drop is 
then set up by, e.g., using a rotating shaft Q] or applying 
an acoustic torque . An elegant way to spin up drops 
made from polarizable fluids is to use rotating electric |^ 
or magnetic fields p^ . 

In the present letter we report theoretical and experi- 
mental investigations of rotating ferrofluid drops and in 
particular perform the first quantitative study of a tran- 
sition from a spheroidal to a non-axissymmetric equi- 
librium shape in this system. Moreover, providing an 
approximate solution of the hydrodynamic flow problem 
inside and outside the drop we are able to analyze the 
rotational motion of the drop shape and to separate it 
from the internal hydrodynamic flow. 

Ferrofluids are suspensions of ferromagnetic nano- 
particles in suitable carrier liquids combining the hydro- 
dynamic behaviour of a Newtonian fluid with the mag- 
netic properties of a super-paramagnet [0. A rotating 
external field induces a rotational motion of the nano- 



particles which due to their viscous coupling to the sur- 
rounding liquid transfer the angular momentum to the 
whole drop. The system has been studied previously 
by Bacri et al. using microdrops with a typical radius 
of 10 /xm, very small surface tension and a viscosity 
much larger than that of the outer fluid j|, Al- 
though an instability of the axissymmetric shape was 
found experimentally and theoretically the emerg- 
ing non-axissymmetric configurations could not be stud- 
ied quantitatively since for the given parameter values 
very irregular shapes arise. In our experiments we used 
much larger drops of radius R = 2.75mm of a kerosine- 
based ferrofluid with magnetic susceptibility /i^ — 15.3 
and dynamic viscosity rji = 0.019 Pas immersed in 3- 
brome-l,2-propandiol with dynamic viscosity 772 = 0.058 
Pas. The interface tension is cr = 2.8 ■ 10"'^ N/m, the 
frequency of the magnetic fleld 560 Hz. For these pa- 
rameter combinations we flnd well-deflned transitions to 
three-axes ellipsoids which can be analyzed in detail. 

From the parameter values given above we find for the 
typical time of shape relaxations of the drop = 0.1s. 
We are thus considering the case of a fast rotating field, 
LOTs ^ 1, in which the drop assumes an oblate shape 
with the short axis perpendicular to the field plane. Its 
form is very near to an ellipsoid and we will use this 
approximation throughout our theoretical analysis. The 
shape will be specified by the ratios = a/6 and Cc = a/c 
between the seminaxes of the ellipsoid where as usual 
a > 6 > c is assumed. 

Estimating the orders of magnitude of the energies due 
to magnetization, surface tension, viscosity and intertia, 
respectively, we find that for the experimentally relevant 
parameter values the shape is almost completely deter- 
mined by the balance between surface and magnetic en- 
ergy where the latter is conveniently averaged over one 
period of the field rotation. The calculation of the mag- 
netic energy is simplified by the fact that an ellipsoid 
in a homogenous external field builds up a homogenous 
magnetization. Using a coordinate system in which the 
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external field is of the form G = (Gcos(wt), Gs'm{u!t),0) 
this magnetization is determined by the relations 



UyMy (1) 



between the components of the external field G, the in- 
ternal field H and the magnetization M |Q. Here 
and Uy denote the demagnetization factors along the x 
and y axis, respectively, which are known functions of 
and Cc 0. 

The calculation of the magnetic energy is most easily 
accomplished assuming a linear magnetization law, M = 
xH, with the susceptibility x = /i^ — 1- We then find for 
the magnetic energy Em = — ^o^M • G/2 using (0) 
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with V denoting the volume of the drop. Averaging over 
one period of the rotation [T5| , using the well-known ex- 
pression for the surface of a three-axes ellipsoid, and ob- 
serving volume conservation for the drop we obtain the 
following result for the sum E of magnetic and surface 
energies 
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X^x 1 + X^y 

{F{m,K) + (4 - l)E{m,K)) 



Here F and E are elliptic integrals of first and sec- 
ond kind, respectively, m = — l/cc, k = 

V i^c ~ ^b)/(^c ~ l)i B = ^oG^R/a is the magnetic 
Bond number measuring the strength of the external 
field. Minimizing this expression numerically in the ge- 
ometry parameters £{, and Ec the dependence of the shape 
of the drop can be determined for varying external field 
strength G. Results for the parameter values given above 
are shown in figs.|l| and ^ together with our experimental 
findings. 

The main result is a transition from a spheroid charac- 
terized by efc = 1 to a pronounced non-axissymmetric 
form of a three-axes ellipsoid for intermediate values 
of the magnetic field strength. This transition takes 
place only if the magnetic permeability is large enough, 
/ir > 5.1. It occurs via supercritical bifurcations up to 
fj,j. = 11.6 and through subcritical transitions with hys- 
teresis for still higher values of fir including our experi- 
mental value. 

From the figures it is seen that the linear theory is in 
qualitative agreement with the experiment. It also yields 
quantitatively good results for the values of the Bond 
number at which the transitions to three-axes ellipsoids 
occur. The ratios between the semiaxes, however, are 
overestimated, with the discrepancy increasing with the 
field strength. In fact the magnetic field corresponding to 




FIG. 1: Ratio et = a/b between the two largest semiaxes of 
a rotating ferrofluid drop as a function of the magnetic Bond 
number B. Squares are experimental values with filled sym- 
bols corresponding to increasing, empty symbols to decreas- 
ing field strength, respectively. The curves show from top to 
bottom the results for a linear magnetization law, M — xH, 
for the Langevin M{H) and for the dynamic curve M{H) as 
determined from an independent experiment. Full lines corre- 
spond to stable configurations, dotted lines to unstable ones. 
There are no fit parameters. 
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FIG. 2: Same as fig.|l] for the ratio ec = a/c between the 
largest and the smallest semiaxes of the rotating drop. 



a Bond number of B = 80 is of the order of the saturation 
magnetization Moo — 80kA/m and hence deviations from 
the linear law M = xH become important. 

In order to account for these deviations from a linear- 
ity we have extended the theory to general magnetization 
curves M{H). The field inside the ellipsoid is still homo- 
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geneous and using (||) we get 
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from which for given M{H) the time dependence of the 
internal field H can be determined numerically. The 
magnetic energy can now be written in the form 
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For a linear magnetization law the last two terms cancel 
and we find back (^. Expression is numerically av- 
eraged over one period of the field rotation and the total 
energy is minimized with respect to eb and ec- 

A standard choice for the magnetization curve of a fer- 
rofluid is the Langevin law M — A/oo(coth(^) — 1/0 with 
^ = 3xH/Maa- The corresponding results are included in 
figs.|l| and||. It is clearly seen that although the values for 
eb and ec are reduced in comparison with the linear the- 
ory the corrections are rather small. The physical reason 
for the remaining discrepancy lies in relaxation effects of 
the magnetization which for the comparetively high fre- 
quency of the external field start to play a role. We have 
therefore determined the magnetization curve M{H) of 
our ferrofluid for an oscillating magnetic field from an in- 
dependent experiment. The lowest curves in figs. I and I 
result from the use of a spline interpolation of these data 
in our numerical code. Whereas the experimental results 
for Ch are now underestimated for large values of B the 
correspondence between theory and experiment for is 
rather satisfactory. 

We now turn to the analysis of the drop motion. In 
the experiments we observe a slow rotation of the non- 
axissymmetric drop shape with angular velocity = 
O.lHz in the direction of the field rotation. In order to 
discuss this motion theoretically we have to solve the 
hydrodynamic flow problem inside and outside the drop 
taking into account the appropriate boundary conditions 
for normal and tangential stresses at the surface. In its 
general form this is a formidable free boundary value 
problem. We hence employ some reasonable approxima- 
tions consistent with our assumption of an elliptical drop 
shape. For the experimentally relevant sizes and viscosi- 
ties the inertial terms in the hydrodynamic equations are 
small and we may use the Stokes approximation for the 
determination of the flow fields. It is convenient to use 
the coordinate system in which the shape of the drop is 
at rest. In this frame we assume that the motion inside 
the drop is of uniform vorticity C and hence use for the 
internal and external flow fields the ansatzes 

V'" = (-C^, C-, 0), V- = V- + {uj, vj, wj). (5) 
a 



From the continuity of the velocities at the surface of the 
drop we infer that {uj,vj,wj) describes the flow field 
outside a rigid elliptical particle at rest in a viscous fluid 
with asymptotic velocity limr^ooiuj,vj,wj) = {{^ + 
(^a/b)y, —{Q + (^b/a)x,0), a problem solved by Jeffery 
many years ago [|l^ . 

To complete the solution in the discussed approxima- 
tion we need two equations to fix the so far undetermined 
parameters and ^. As first equation we use the bal- 
ance between the viscous torque experienced by the drop 
in its rotational motion and the magnetic torque, again 
averaged over one period of the driving. The second one 
results from the stationarity of the energy balance stat- 
ing that the average work done by the external field per 
unit time must be equal to the energy dissipated per unit 
time in the viscous flows. This latter reqirement can 
be rewritten as a continuity condition for the tangential 
stress averaged over the surface of the drop |l^ . 

The expression for the viscous torque builds on Jef- 
fery's solution and can be obtained in the same way as 
for a slowly rotating field, luTs ^ 1, ||2^. For the deter- 
mination of the magnetic torque L = ^qV (M x G) it is 
crucial to take into account the finite relaxation time of 
the magnetization since a non-zero averaged torque re- 
sults from the phase lag between M and H. The explicit 
calculation is again feasible analytically for a linear mag- 
netization curve M = xH where the susceptibility is now 
complex, X = Xi ~ *X2- We will assume X2 ^ Xi consis- 
tent with the experiment where xi — 14.3 and X2 — 2.1. 
Using (|l|) and uj — Q = lu we get for the z-component of 
the averaged magnetic torque 

— _ t^oVG' ( 1 , 1 ^ 
- \{l + Xn.? + [l + xny?) ■ 

In the integrated balance of tangential stresses the vis- 
cous contributions can be obtained from Jeffery's solution 
|l7| , [l^ . The tangential part of the Maxwell stress ten- 
sor of the magnetic field is automatically continuous at 
the interface due to the boundary conditions obeyed by 
the fields |l|. However, in the present non-equilibrium 
situation there is a magnetic contribution due to the ad- 
ditional term (MiHk—MkHi) /2 in the stress tensor intro- 
duced by Shliomis |^l| . Taking into account the various 
contributions two equations for and C can be derived. 
The final expressions are rather long and will be pub- 
lished elsewhere p^ . 

Results for the rotation frequency H. of the drop are 
shown in fig]^ together with the corresponding exper- 
imental findings. The linear theory overestimates the 
magnetic torque and therefore also the frequency of ro- 
tation. To investigate the effects of deviations from the 
linear magnetization curve we have calculated the mag- 
netic torque using the numerical solution of the magne- 
tization relaxation equation [ p3[ 

— = ^(H • (M - Mo))H - -l^H X (M X H). 
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For larger values of the Bond number the peaks get 
smaller and more numerous and the approximation of 
the shape by an ellipsoid becomes again more accurate 
resulting in a reasonable agreement between theory and 
experiment. 
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FIG. 3; Rotation frequency of the drop as function of the 
magnetic Bond number. Squares are experimental values, 
lines are theoretical results using a linear magnetization law 
(dashed) and the magnetization curve M{H) as determined 
in an independent experiment (full), respectively. The up- 
per part describes the rotation of the three-axes ellipsoid, the 
lower one the hard body rotation of a disk. 



The static magnetization Mo(-ff) and the field depen- 
dence of the relaxation times r,, and t± are determined 
from independent experiments, the geometry of the drop 
is taken from figs.|^ and ^. Although the magnetic torque 
is substantially smaller than in the linear case the re- 
sults for the rotation period of the non-axissymmetric 
drop are almost the same. This is due to the reduced 
viscous torque resulting from smaller values of eh and ec 
in the nonlinear theory. The remaining differences with 
the experiment may be due to the polydispersity of the 
fluid requiring a whole spectrum of relaxation times and 
a small ellipticity of the external field. These questions 
will be dealt with in detail elsewhere p2) . 

For the axissymmetric state, ef, = 1, shown in the 
lower part of fig.§ the situation is somewhat simpler. In 
this case the drop forms a flat disk rotating like a hard 
body. For large values of B this disk developes peaks 
around its perimenter due to the normal field instability 
^] making an experimental determination of the 
rotation frequency easy. As can be seen from fig.^ the 
linear theory again strongly overestimates the rotation 
frequency whereas the non-linear theory yields substan- 
tially smaller results. For Bond numbers between 80 and 
100 the disk shows few, rather large peaks resulting in an 
increased viscous torque and therefore a slower rotation. 
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